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Tank  Physics 


This  report  discusses  the  quality  of  SGEMP  simulation  which  would  exist  in 
the  Satellite  X-ray  Test  Facility  if  it  were  constructed  in  accordance  with 
the  proposed  strawman  design. 

It  was  found  that  the  simulation  quality  at  some  locations  on  a  test  satellite 
is  poor,  suggesting  that  a  modification  of  the  strawman  design  may  be 
necessary. 
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SECTION  1 
INTRODUCTION 


It  is  the  purpose  of  this  report  to  present  the  results  of  a 
recent  investigation  into  the  quality  of  SCHLMP  simulation  which  can  be 
expected  in  the  proposed  Satellite  X-rav  Test  facility.  The  strawman  design 
which  will  he  criticized  in  this  report  consists  of  a  spherical  vacuum  tank 
30  m  in  diameter  with  a  concentric,  transparent,  non-emitting  damper  membrane 
whose  diameter  is  80  percent  of  the  tank's  diameter.1  The  X-ray  illumination 
is  caused  by  both  a  compact  exploding  wire  and  an  extended  bremsst rahlung 
source.  In  this  study,  we  wish  to  determine  how  closely  the  SGPMP  response 
of  a  test  body  placed  in  this  simulator  corresponds  with  the  response  of  the 
same  boils  in  free  space. 


In  order  to  examine  the  simulation  quantitatively,  we  compare 
tile  results  of  code  calculations  of  the  StlfMP  response  of  a  test  body  placed 
in  the  tank  to  that  of  the  same  body  in  free  space.  Several  potentially 
important  effects  are  introduced  by  the  tank.  first,  electromagnetic  energy 
cannot  he  radiated  into  free  space  but  is  instead  reflected  from  the  tank 
walls  and  returned  to  the  test  object.  These  reflections  cause  the 
deposition  of  additional  electromagnetic  energy  in  the  satellite  possibly 
causing  severe  overtesting  and  even  satellite  damage.  In  an  effort  to 
overcome  this  difficulty  the  resistive  membrane  is  placed  inside  the  tank. 

An  additional  complication  is  caused  by  the  emission  of  electrons  from  the 
tank  walls.  These  particles  are  a  source  of  unwanted  electromagnetic  fields 
and  can  also  charge  the  test  object. 
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In  order  to  isolate  each  of  these  effects  we  have  performed  calcu¬ 
lations  in  which  no  electron  emission  occurs  from  the  walls  or  from  the 
membrane.  This  calculation  simulates  perfect  suppression  of  wall  and  membrane 
emission.  It  demonstrates  the  effects  of  reflection  from  the  tank  walls  and 
absorption  of  energy  by  the  resistive  membrane.  Ke  then  examine  i he  quality 
of  simulation  provided  by  the  strawman  design  in  which  electron  emission 
from  the  tank  is  included,  finally,  we  examine  the  extreme  case  of  an 
emitting  membrane  and  emitting  wall.  Some  damper  designs  have  proposed  the 
use  of  sheets  of  conducting  low  Z  material  rather  than  ohmic  meshes  and 
in  this  calculation  we  study  the  effects  of  membrane  emission.  Although  this 
calculation  may  not  prove  to  he  relevant  to  the  SXTF  design,  it  is  reported 
here  since  it  may  he  of  interest  in  connection  with  other  SGEMP  experiments. 

It  is  the  aim  of  these  studies  to  determine  whether  the  passive 
suppression  provided  by  the  low  Z  emission  from  the  walls  of  the  tank  is 
adequate  or  whether  active  suppression  is  required. 

Since  the  SXTF  may,  at  some  future  date,  he  used  to  test  satellites 
whose  dimensions  are  comparable  to  the  dimensions  of  the  tank,  we  wish  to 
determine  to  what  extent  the  SGI- Ml’  is  affected  when  parts  of  the  satellite 
are  in  close  proximity  to  the  tank  walls  and  to  the  membrane. 

In  Section  2  we  describe  our  method  of  calculation  and  discuss  the 
various  approximations  and  assumptions  made.  A  detailed  presentation  of  the 
results  of  our  calculations  is  provided  in  Section  3.  Section  4  is  devoted 
to  a  discussion  of  the  significance  of  our  results  and  in  Section  5  we  make 
our  recommendations  for  modifications  of  the  strawman  design. 
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SECTION  2 

METHOD  OF  CALCULATION 


Our  calculations  were  performed  using  the  two-dimensional 
eyl  indrical ly  symmetric  particle  pushing  code  SF.MP,  a  detailed  account  of 
which  is  provided  in  Reference  2.  Several  important  modifications  to  this 
code  were  required  for  our  tank  simulations  and  we  shall  discuss  these 
modifications  in  this  section. 

2.1  CODE  DESCRIPTION  OF  THE  SXTF 

The  tank  simulator  is  sketched  in  Figure  1.  Basically,  it  con¬ 
sists  of  a  spherical  vacuum  tank  (stainless  steel )  15  m  in  radius  from  which 
a  segment  has  been  removed  in  order  to  provide  a  suitable  location  for  the 
radiation  sources.  The  total  length  of  the  assembly  is  thus  24  m.  In  SUMP, 
this  is  represented  by  a  grid  of  155  s  85  cells  each  18  cm  in  the  radial 
and  axial  directions.  A  resistive  membrane,  spherical  in  shape  and  of  radius 
12  meters  is  positioned  concentric  with  the  sphere  and  the  satellite  was 
placed  inside  the  membrane. 

2.2  MESH  AND  BOUNDARY  CONDITIONS 

As  with  any  finite  difference  scheme,  our  mesh  must  be  small 
enough  to  provide  an  adequate  description  of  the  phenomenon  of  interest,  in 
our  case  the  motion  of  the  electrons  and  the  fields  so  produced,  yet  large 
enough  so  that  calculation  does  not  become  excessively  time  consuming  or 

require  unduly  large  amounts  of  core.  The  mesh  size  determines  the  time 

step  because  of  the  requirement  that  Courant's  condition  be  satisfied. 
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Tank 


Figure  1.  Schematic  drawing  of  SXTF.  The  membrane  is  3  m 
from  the  tank  walls. 
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IV  e  therefore  require  as  large  a  mesh  as  possible  consistent  with  accurate 
representation  of  the  motion  of  charge. 


The  original  STMI'  mesh'’  used  three  interlocking  meshes  with  the 
bod\  drawn  in  the  smallest  mesh  and  the  outer  boundary  in  the  largest  mesh. 
Our  free  space  calculations  used  the  three  meshes  and  the  transmitting  outer 
boundary  condition.  In  our  tank  calculations,  however,  we  have  a  reflecting 
and  emitting  outer  boundary,  the  walls  of  the  tank.  Accurate  representation 
of  the  emission  of  particles  from  the  tank  walls  demands  that  the  mesh  spac¬ 
ing  at  the  outer  boundary  be  smaller  than  the  Debye  length  in  the  charge 
cloud.  Satisfaction  of  this  condition  precludes  the  use  of  the  three  mesh 
system  and  a  single  uniform  mesh  was  used. 

A  sphere  was  drawn  inside  this  mesh.  The  usual  boundary  conditions 
were  applied  at  the  surface  of  the  sphere.  The  walls  of  the  tank  are  assumed 
to  be  infinitely  conducting  so  that  the  tangential  component  of  the  electric 
field  must  vanish  at  its  surface.  The  reflecting  boundary  condition  is  also 
applied  here. 


The  mesh  spacing  in  these  calculations  was  uniform  with  Ar  =  At  = 

IS  cm.  Under  conditions  of  extreme  space-charge- 1 iini ting,  this  mesh  is 
inadequate  since  the  Debye  length  is  smaller  than  the  mesh  spacing  (X^  ~ 

1  cm).  This  condition  prevails  at  the  emitting  surface  of  the  satellite 

but  not  at  the  tank  walls  | \  -  JO  cm).  At  the  emitting  surface  of  the 

satellite,  particle  motion  cannot  be  described  in  this  mesh  and  we  must  use 
an  alternative  description  of  the  currents  a t  that  position.  In  a  space- 
charge- 1  i  mi  t  ed  situation,  a  dipole  layer  forms  above  an  emitting  surface, 
ihe  current  density  may  be  described  in  terms  of  the  time  derivative  of  the 
dipole  moment  per  unit  area.  A  P  driver  may  he  calculated  for  a  given 

X-ray  spectrum  and  time  history  for  a  given  material  using  the  one-dimensional 

particle  pushing  code  SUM, ID,  details  of  which  may  he  found  in  Reference  5. 
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In  the  SXTf'  strawman  design,  two  sources  of  radiation  are  used  —  a 
compact  exploding  wire  source  and  a  more  extended  hremsstrahlung  source. 

Since  the  hremsstrahlung  spectrum  is  considerably  harder  than  the  exploding 
wire  spectrum  (2  keV  hlackbody),  fewer  electrons  are  produced  so  that  the 
Debye  length  in  the  charge  cloud  produced  by  the  hremsstrahlung  is  larger 
than  the  cell  size.  l-'or  this  reason,  a  dipole  driver  description  of  the 
hremsstrahlung  excitation  is  unnecessary,  although  it  is  required  to  simulate 
the  boundary  layer  associated  with  the  exploding  wire  photon  source. 

2.3  OTHER  CODE  MODIFICATIONS 

In  the  SEMI'  code,  particles  inside  a  particular  cell  are  moved 
using  interpolated  values  of  the  electric  and  magnetic  fields.  The  interpola¬ 
tion  is  linear  using  the  field  values  at  the  corners  of  the  cell.  If  one  has 
a  conductor  (e.g.,  the  damping  membrane)  embedded  in  the  mesh,  the  fields  are 
discontinuous  at  its  surface.  Under  these  circumstances,  a  reliable  estimate 
of  the  field  at  a  particular  point  can  no  longer  be  obtained  by  interpolation 
and  we  use  the  nearest  grid  point  method.  In  this  method,  particles  are 
moved  using  the  values  of  the  fields  at  the  grid  point  closest  to  the  particle, 
on  the  same  side  of  the  conductor.  One  thus  avoids  the  difficulty  of  inter¬ 
polating  across  a  d i scont i nu i ty . 

In  our  calculations,  we  have  represented  a  spherical  surface  in  a 
cylindrical  mesh.  The  cylindrical  mesh  was  chosen  as  a  matter  of  convenience 
since  the  satellite  is  more  easily  represented  in  a  cylindrical  mesh.  This 
representation  of  a  spherical  surface  can  lead  to  inaccuracies  in  the  rate 
of  propagation  of  fields  if  the  wall  of  the  tank  [Kisses  diagonally  through 
a  cell.  Consider  the  cell  shown  in  figure  2  through  which  the  tank  wall 
passes  diagonally.  A  signal  propagating  from  A  to  11  reaches  B  in  a  time 
VT'.A/c,  in  the  code  calculations,  such  a  signal  would  be  required  to  travel 
through  C,  reaching  B  in  a  time  2\r/c. 
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In  order  to  correct  this  inaccuracy,  we  proceed  as  follows. 
Faraday' s  equation  may  be  written 


cel  1 


In  any  ordinary  cell  of  the  mesh,  we  can  write  the  integral  on  the  left  as 

2  I A  r  A :: )  li  , 

but  when  the  tank  wall  passes  diagonally  through  a  cell,  the  integral  is 
over  a  half  cell  and  we  have 


Using  Stoke 's 


Removing  the 


ADB 


theorem,  we  have 

integration  signs,  we  obtain 


•  da  . 


Us  i  ng 


Amperes  law  may  tie  written 

rxv:  rxv; 

Stoke's  theorem  and  removing  the 

V  '<  B  =  -In. J  +  --  . 

~  ~  "  c  at 


integration  sign  we  find  that 


Since  the  integration  is  over  the  cell  TXYZ,  Ampere's  equation  is  unchanged 


Combining  Ampere's  and  Faraday's  equations,  we  have 


V  *  V  x  F 
In  vacuum,  we  have 


v2f  =  J_ 


13,..  1  3E, 

+  r  iff)  ■ 


2c  3t 


3t“ 

This  is  a  wave  equation  with  effective  propagation  velocity  V2c.  The 
introduction  of  the  factor  of  2  in  Faraday's  equation  corrects  the  propa¬ 
gation  velocity  in  those  cells  through  which  the  tank  wall  passes  diagonally. 
In  all  other  cells  this  correction  is  unnecessary. 


2.4 


RETARDATION  OF  THE  RADIATION  PULSE 


Since  the  duration  of  the  X-ray  pulse  from  both  the  exploding 
wire  source  and  the  bremsstrahlung  source  is  shorter  than  the  light  travel 
time  across  the  tank,  care  must  be  taken  to  retard  the  emission  from  a  point 
by  the  appropriate  light  travel  time.  The  zero  point  of  time  occurs  when 
the  first  emission  point  (either  on  the  satellite  or  on  the  wall  of  the 
tank)  starts  to  emit.  The  first  emission  point  is  the  emission  point  closest 
to  the  radiation  source.  Let  its  2  coordinate  be  denoted  by  Z^.  A  point  at 
2^  will  not  begin  to  emit  until  the  pulse  has  reached  it,  i.e.,  at  time 
Z^-Z()/c.  Every  emission  point  emits  with  some  pulse  time  history.  However, 
the  pulse  at  each  point  is  displaced  by  the  amount  (Z-Z^)/c. 

The  X-ray  pulse  is  thus  treated  as  a  plane  wave  propagating  in 
the  Z  direction;  the  time  of  emission  from  a  point  is  determined  by  its  Z 
coordinate.  The  plane  wave  approximation  is  appropriate  when  we  are  con¬ 
sidering  the  extended  bremsstrahlung  source  but  may  be  slightly  in  error 
when  we  are  considering  the  more  compact  exploding  wire  source  which  is 
better  represented  by  a  spherical  wave. 
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2.5 


DAMPING  MEMBRANE 


The  membrane  serves  two  purposes.  In  all  SXTF  ealculat ions ,  the 
membrane  extracts  energy  from  the  wave  reflected  from  the  tank  walls.  Its 
purpose  is  to  diminish  the  effect  of  reflection  so  that  the  tank  more  closely 
simulates  free  space.  in  some  calculations,  we  examined  a  non-transparent 
membrane  which  emits  electrons  and  through  which  particles  cannot  pass.  In 
these  calculations,  the  membrane  has  the  additional  functions  of  electron 
emitter  and  absorber. 

The  resistive  membrane  is  simulated  by  assigning  a  resistance  to 
those  cells  which  correspond  to  the  position  of  the  membrane.  It  has  been 
shown4  that  the  optimal  membrane  impedance  is  I  =  Z^/2  where  is  the 
impedance  of  free  space  .378  Since  the  behavior  of  the  fields  is  not  very 
sensitive  to  the  value  of  the  impedance,  we  chose  a  value  of  200  ft/square. 

We  then  reduced  the  electric  field  components  by  a  factor  of  1/(1  +  Z^cAt/ 
Z.\)  where  c  is  the  speed  of  light  and  At  and  A  are  the  time  step  and  spatial 
size  respectively.  ibis  procedure  removes  energy  from  the  tank  in  order  to 
produce  a  better  simulation  of  free  space. 

2.6  THE  TEST  SATELLITE 

The  satellite  chosen  for  our  simulation  test  was  a  FLTSATCOM  like 
the  model,  a  rough  sketch  of  which  is  shown  in  Figure  3 a .  This  choice  was 
made  since  it  is  of  interest  to  examine  the  response  of  a  resonant  body. 

It  may  stress  the  simulation  fidelity  requirement  to  a  greater  extent  than 
a  non-resonant  body  because  it  has  characteristic  frequencies  which  persist 
in  time. 

Our  satellite  model  shown  in  Figure  3b  consists  of  three  cylinders 
each  of  length  2  m  and  radius  1  m  connected  by  thin  struts  2.S  cm  in 
diameter.  Since  the  struts  cannot  be  properly  described  in  this  mesh,  a 
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special  treatment  was  necessary.  Our  treatment  of  struts  is  outlined  in 
Appendix  A. 

2.7  PARAMETERS  OF  THE  CALCULATIONS 

We  now  give  the  numerical  values  of  the  various  quantities  used 
in  our  calculations. 

All  calculations  presented  in  Section  3  assume  an  X-ray  fluence  of 
-4  ■> 

5  x  10  cal/cm“/sec.  The  X-ray  flux  has  two  distinct  components  —  a  brems- 
strahlung  component  and  an  exploding  wire  component.  The  bremsstrahlung 
component  has  an  essentially  flat  spectrum  as  shown  in  Figure  4a  and  a 
triangular  time  history  which  peaks  in  15  ns  and  falls  slowly  to  zero  at 
70  ns  (Figure  4b).  The  exploding  wire  spectrum  is  that  of  a  2  kcV  black 
body  and  its  time  history  is  somewhat  complicated  (Figure  4c).  A  peak  of  the 
exploding  wire  radiation  occur s  after  only  7.5  ns  and  is  followed  by  a  rapid 
decline.  This  initial  peak  is  followed  by  several  secondary  peaks.  Because 
of  Dehye  length  considerations  discussed  previously,  exploding  wire  radiation 
is  represented  at  the  surface  of  the  satellite  by  the  l1  driver  shown  in 
Figure  4d . 


The  energy  spectrum  of  electrons  emitted  from  the  satellite  was 
taken  to  be  the  backward  emitted  spectrum  for  aluminum.  Emission  from  the 
tank  assumes  the  forward  emitted  spectrum  of  carbon.  Use  of  the  forward 
spectrum  is  not  strictly  correct  but  comparison  of  the  forward  and  backward 
spectra  of  carbon  shows  that  the  error  introduced  is  small.  Although  the 
photoelectric  yield  of  curb on  is  smaller  than  that  of  aluminum,  the  walls 
of  the  tank  have  a  much  larger  area  than  the  satellite  and  consequently 
emit  approximately  an  order  of  magnitude  more  electrons. 
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Figure  4.  Spectra  and  time  histories  of  X-ray  drivers. 
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SECTION  3 

RESULTS  OF  THE  CALCULATIONS 


The  present  section  is  the  first  of  two  devoted  to  a  detailed 
presentation  and  discussion  of  the  results  of  our  calculations.  A  brief 
summary  of  the  initial  assumptions  of  the  various  runs  is  given  in  Table  1. 

3.1  SATELLITE  IN  FREE  SPACE 

Since  in  a  perfect  simulator,  the  SGHMP  response  of  a  body  inside 
a  test  tank  should  he  the  same  as  that  of  the  same  body  in  free  space,  it  is 
useful  for  comparison  purposes  to  calculate  the  free  space  response  of  our 
test  body. 


In  the  free  space  run,  the  outer  boundary  condition  is  such  as  to 
allow  essentially  all  of  the  electromagnetic  energy  reaching  the  outer 
boundary  to  be  radiated.  The  treatment  of  the  outer  boundary  condition  in¬ 
troduces  a  fractional  error  of  approximately  \/2'tR  (see  Reference  5),  where 
\  denotes  the  longest  wavelength  of  electromagnetic  radiation  of  interest 
(which  is  of  the  order  of  the  length  of  the  satellite-)  and  R  is  the  radial 
distance  to  the  outer  boundary,  for  accuracy,  it  is  required  that  the 
fraction  X/d'R  be  small  compared  to  unity.  This  condition  requires  that  the 
distance  to  the  outer  boundary  be  much  larger  than  the  dimensions  of  the 
test  body.  Computer  storage  limitations  make  it  difficult  to  satisfy  this 
requirement  while  providing  a  sufficient  number  of  cells  for  an  accurate 
description  of  the  satellite.  To  satisfy  this  requirement  while  keeping 
computational  demands  at  a  reasonable  level,  we  made  use  of  the  three 
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interlocking  meshes  of  increasing  sice  of  the  original  SUMP  code  rather 
than  using  the  single  mesh. 

In  a  perfect  SGEMP  simulator,  the  fields  should  be  identical  to 
those  calculated  for  the  frec-space  ease  being  free  from  reflected  signal 
and  oscillations  resulting  from  the  excitation  of  the  natural  modes  of  the 
tank . 


An  additional  effect  is  introduced  by  the  /iresence  of  the  tank. 
Although  the  walls  of  the  tank  are  coated  with  carbon,  a  material  of  low  Z 
and  low  photoelectric  yield,  the  tank  has  a  much  larger  surface  area  than  the 
satellite  and  emits  considerably  more  electrons.  The  irradiation  of  the 
tank  and  satellite  by  X  rays  produces  a  larger  quantity  of  electromagnetic 
energy  than  was  produced  by  the  irradiation  of  the  satellite  in  free  space 
and  it  is  of  interest  to  examine  the  effects  of  this  additional  energy. 

3.2  OVERVIEW  OF  THE  CALCULATIONS 

In  order  to  examine  quantitatively  the  effects  of  reflection  and 
photoelectric  emission  from  the  tank,  three  other  cases  have  been  calculated. 

Case  A— Perfect  Suppression 

In  this  case,  no  emission  takes  place  cither  from  the  walls  of 
the  satellite  or  from  the  membrane.  The  only  emitting  surface  in  the  problem 
is  the  front  surface  of  the  satellite.  Wo  call  this  case  perfect  suppression 
because  the  electron  emission  from  the  tank  walls  are  completely  suppressed. 

The  assumption  of  a  non-emitting  tank  is  not  physically  justifiable 
but  the  calculation  is  intended  to  model  perfect  suppression  of  electron 
emission  from  the  tank  walls  and  isolate  the  effects  of  reflection  from  the 
tank . 
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Case  B— Strawman 


In  this  case,  we  model  the  strawman  design  for  the  SXTF  described 
in  Reference  1.  The  satellite  is  placed  inside  a  vacuum  tank  whose  walls 
are  coated  with  carbon.  Emission  takes  place  from  the  walls  with  the  photo¬ 
electric  yield  appropriate  to  carbon.  The  damper  membrane  is  described  in 
Reference  1  as  an  ohmic  mesh  made  up  of  carbon  coated  fiberglass  0.5  cm  in 
diameter  made  into  a  mesh  with  10  cm  spacing.  In  the  calculations  we  assign 
the  appropriate  resistance  to  the  membrane,  200  0/square,  and  assume  that  it 
is  non-emitting  and  transparent  in  the  sense  that  electrons  pass  freely  through 
it. 


This  calculation  is  intended  to  be  physically  realistic  and  is 
designed  to  isolate  the  potentially  important  effects  of  electron  emission 
from  the  walls. 

Case  C— Emitting  Walls  and  Membrane 

Since  other  designs  of  damper  membranes  have  been  considered  for 
the  SXTE,  we  examine  the  case  where  a  non-transparent,  emitting  carbon  membrane 
is  used.  In  this  case,  electrons  are  not  allowed  to  pass  freely  through  the 
membrane  but  are  absorbed  upon  reaching  it.  The  membrane  also  emits  from 
both  sides  with  the  photoelectric  yield  of  carbon. 

Cases  A-C  were  calculated  with  the  satellite  centered  and  again 
with  the  satellite  close  to  the  emitting  membrane.  The  latter  (shifted) 
configuration  was  designed  to  test  the  situation  in  which  a  part  of  the 
satellite  was  close  to  an  emitting  surface. 
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3.3 


OVERVIEW  OF  THE  RESULTS 


A  useful  measure  of  the  .Joule  heating  of  electronic  circuitry  due 
to  an  induced  electromotive  force  is  the  quantity  /B^ftjdt.  Comparison  of 

this  quantity  for  various  cases  provides  an  indication  of  the  simulation 

fidelity.  We  are  also  interested  in  the  behavior  of  the  skin  current  and 

normal  components  of  the  electric  fields  at  various  points  on  and  near  the 

surface  of  the  satellite. 

In  this  report  all  electric  fields  are  expressed  in  statvolts 
per  centimeter.  Skin  currents  arc  reported  in  Amperes  per  meter.  The  time 
integral  of  the  square  of  the  magnetic  field  is  always  given  in  arbitrary  units. 

The  fields  are  calculated  and  plotted  for  ten  points  but  here  we 
present  the  results  for  two  representative  points.  Observer  point  3  (see 
figure  3)  is  close  to  the  emitting  surface  of  the  satellite.  Observer 
point  10  is  on  the  surface  farthest  from  the  radiation  source.  The  signal 
at  point  3  is  dominated  by  the  emission  from  the  front  face  of  the  satellite. 

At  point  10,  the  effects  of  reflected  signals  and  electron  emission  from  the 
tank  arc  much  more  readily  seen.  Our  main  results  may  be  summarized  as  follows. 

1  he  simulation  fidelity  is  adequate  at  point  .S  satisfying  the  fidelity 
criteria  described  in  Relerence  (>.  Departures  from  the  free  space  response 
occur  in  both  the  strawman  and  emitting  membrane  cases.  The  electric  field 
at  point  10  at  late  time  is  dominated  by  a  tank  induced  pulse  which  is  stronger 
than  the  initial  pulse.  Considerably  more  energy  is  deposited  at  observer 
point  10  in  the  strawman  and  emitting  membrane  cases  than  in  the  free  space 
case.  At  observer  10,  the  unshifted  satellite  strawman  case  fails  to  satisfy 
the  criteria  for  adequate  simulation  fidelity. 

In  the  case  where  the  satellite  is  positioned  off  center,  the 
departure  from  the  tree- space  response  is  more  severe.  Simulation  fidelity 
is  still  adequate  at  observer  point  3  but  at  point  10  dramatic  differences 
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appear.  The  behavior  of  the  electric  field  at  point  10  is  characterized  by 
a  large  negative  spike  occurring  at  the  time  of  reflection.  The  spike  is 
an  order  of  magnitude  greater  than  the  initial  spike  in  the  free  space  case. 

The  integral  of  B“  is  greater  than  in  the  free  space  case  by  a  factor  of  200. 
The  simulation  fidelity  at  point  10  in  the  shifted  satellite  case  must  be 
cons i de  red  i nadcquat e . 

We  do  not  consider  that  an  overtest  has  occurred  since  the  fields 
at  point  10  in  the  strawman  case  arc  comparable  to  those  at  point  3  for  the 
free  space  case.  However,  the  fields  at  point  10  are  not  faithfully  reproduced 
by  the  simulator. 

3.4  SIMULATION  QUALITY  FOR  CENTERED  SATELLITE  CASE 

We  now  present  the  results  of  our  test  simulations  and  discuss 
in  detail  the  quality  of  simulation  paying  particular  attention  to  the 
causes  of  the  departures  of  the  SGF.MP  response  from  the  free-spaee  response. 

We  first  examine  those  test  simulations  in  which  the  satellite  was  centered 
in  the  tank. 

In  Figure  S,  we  compare  the  integral  of  B~  for  all  four  cases  at 
observer  point  3.  The  general  behavior  of  the  integral  is  similar  in  all 
cases  rising  rapidly  from  zero  and  flattening  off  at  late  time.  Slight 
differences  appear  in  the  perfect  suppression  and  strawman  cases  and  slightly 
higher  l~  20  percent)  values  of  the  integral  are  attained.  This  indicates 
that  reflection  from  the  tank  walls  causes  the  deposition  of  additional 
energy  at  observer  point  3.  The  differences  between  Figures  5b  and  Sc  tire 
attributable  to  electron  emission  from  the  walls  of  the  tank.  Reflection 
from  the  tank  walls  produces  a  slight  ~  10  percent  increase  in  the  maximum 
value  of  ,-B“(t)dt  and  electron  emission  produces  a  further  20  percent  increase. 
\s  one  would  expect,  the  early  parts  of  the  curves  ar<'  almost  identical 
with  the  differences  occurring  at  late  time  t  ■>  100  ns.  These  figures 
indicate  that  noticeably  larger  amounts  of  energy  are  deposited  at  observer 
point  3  in  the  tank  simulation  cases  than  in  the  free  space  case,  but  the 
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Free  Space 


b.  Perfect  Suppression 


effects  of  reflection  and  emission  from  the  tank  walls  do  not  appear  to  be 
serious  at  observer  point  3. 


The  integral  is  plotted  in  figure  Sd  for  the  case  in  which  the 
membrane  emits  and  absorbs  electrons.  Comparison  of  Figure  5c  and  5d 
indicates  that  the  absorption  and  emission  of  electrons  by  the  damper  mem¬ 
brane  has  no  major  effect  on  the  quantity  of  energy  deposited  at  this 
observer  point. 

We  now  examine  the  fields  and  currents  produced  in  these  tests. 
Figure  6  shows  the  skin  current  at  observer  3.  Little  difference  is  seen 
between  the  various  cases  at  point  3;  each  shows  a  large  initial  spike, 
followed  by  subsequent  oscillations.  In  the  free-space  case,  the  skin 
current  had  returned  to  zero  by  100  ns  whereas  oscillations  continued  at 
late  time  when  the  satellite  was  placed  inside  the  vacuum  tank.  The  straw- 
man  and  emitting  membrane  cases  reveal  that  the  late-time  oscillations  are 
stronger  when  more  electrons  are  present  (compare  Figure  6b  with  Figures 
6c  and  d) . 


Figure  7  shows  the  behavior  of  the  normal  component  of  the  electric 
field  at  observer  point  3.  The  changes  in  the  normal  component  of  the 
electric  field  at  observer  3  resulting  from  reflection  or  emission  from  the 
tank  are  very  small  at  observer  point  3. 


From  the  results  shown  in  Figures  5  through  7,  we  must  conclude 
that  the  quality  of  simulation  at  observer  point  3  in  the  centered  satellite 
cases  is  very  good.  Because  observer  point  3  is  close  to  the  emitting 
surface,  the  magnitude  of  the  reflected  pulse  is  small  compared  to  the 
initial  pulse.  It  is  for  this  reason  that  little  difference  is  seen  between 
the  tree-space  case  and  the  case  in  which  the  satellite  is  enclosed  in  the 
tank  . 
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A  more  dramatic  demonstration  of  the  consequences  of  reflection 
from  the  walls  of  the  tank  is  observed  at  point  10.  figure  8  shows  the 
behavior  of  /B“dt  as  a  function  of  time  at  observer  point  10  for  the 
centered  satellite  cases.  Here  the  general  behavior  of  the  fields  is  not 
dominated  by  the  satellite  emission  as  it  was  at  point  3.  In  the  plots  shown  in 
figure  8,  one  can  see  clearly  the  signal  reflected  from  the  tank  walls.  At 
this  point,  the  reflected  pulse  is  actually  larger  than  the  initial  pulse 
from  the  satellite.  We  point  out  that  these  curves  do  not  show  the  total 
effect  of  the  reflected  pulse  since  the  calculat ions  were  terminated  at  a 
time  when  the  integral  was  still  increasing  rapidly  (see  figures  8b,  c  and  d) . 

Because  of  its  distance  from  the  emitting  surface,  the  initial 
pulse  deposits  much  less  energy  at  observer  point  10  than  at  observer  point 
.3 (two  orders  of  magnitude  less,  see  figures  5a  and  8a.  In  the  perfect 
suppression  case,  the  maximum  value  of  the  integral  is  an  order  of  magnitude 
higher  than  in  the  f roe-space  case.  This  increase  is  purely  the  result  of 
reflection  from  the  tank  walls,  figure  8c  indicates  that  the  emission  of 
electrons  from  the  tank  walls  produces  a  further  factor  of  3  increase  in 
the  integral  at  late  time  t  >  100  ns.  fmission  and  absorption  of  electrons 
from  the  membrane  has  little  additional  effect  on  the  energy  deposited  at 
observer  point  10  (see  Figures  5d  and  8d) . 

from  these  results,  we  find  that  using  the  strawman  design  may 
cause  differences  of  a  factor  of  30  in  the  integral  of  B“  at  observer 
point  10  (compare  figures  8a  and  8c).  Since  the  integral  provides  a  measure 
of  the  Joule  heating  of  the  internal  circuitry  of  the  satellite,  wc  must 
conclude  that  the  influence  of  the  tank  walls  may  cause  large  increases  in 
the  amount  of  energy  received  at  observer  point  10.  We  may  add  that  our 
estimate  of  a  factor  of  50  increase  may  be  an  underestimate  since  the 
calculations  were  terminated  at  a  time  when  the  value  of  the  integral  was 
increasing  rapidly.  A  satellite  placed  inside  a  test  tank  could  continue 
to  extract  energy  from  the  cavity  causing  the  energy  received  at  certain 
points  to  exceed  its  free  space  value  by  a  large  factor. 
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It  is  not  clear  that  use  of  the  strawman  design  for  the  SXTF 
would  lead  to  overtesting  of  this  type  of  satellite.  Although  point  10 
receives  much  more  energy  in  the  strawman  case  than  it  docs  in  the  free-space 
case,  the  total  amount  of  energy  deposited  at  point  10  is  much  less  than 
that  deposited  at  point  3.  Since  a  satellite  must  he  able  to  withstand 
irradiation  from  any  direction,  this  cannot  be  considered  an  overtest. 
However,  at  the  particular  observer  point  under  discussion,  this  test  fails 
to  meet  the  criteria  for  simulation  fidelity  which  requires  that 


where  the  subscripts  sim  and  thr  refer  to  the  simulation  and  threat  sources 
respectively.6  (We  have  substituted  the  free  space  response  for  the  threat 
source  in  the  above  relation.) 


We  now  consider  the  fields  and  currents  produced  at  observer  point 
10.  Figure  9  shows  the  behavior  of  the  skin  current  at  observer  point  10 
for  all  four  cases.  Here  obvious  differences  are  seen  in  tin-  tank  simulation 
cases  as  compared  to  the  free-space  case.  The  initial  pulse  at  observer  10 
is  much  smaller  (by  a  factor  of  10)  than  the  initial  pulse  at  observer  point 
3.  While  the  initial  pulse  at  point  10  is  of  the  same  sire  in  all  four 
cases,  the  reflection  causes  large  differences  to  occur  at  late  time.  At 
this  point,  the  reflected  pulse  is  much  larger  than  the  initial  pulse  by  a 
factor  of  2  for  the  case  of  perfect  suppression  and  a  factor  of  3  for  the 
strawman  case  and  for  the  case  of  the  emitting  membrane.  Here  the  emission 
of  electrons  has  a  considerable  effect  on  the  amplitude  of  the  reflected 
pu 1 se. 


From  Figure  9  we  see  that  at  late  time  (t  ''  100  ns),  the  magnitude 
of  the  skin  current  is  one  order  of  magnitude  larger  in  the  strawman  case 
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than  in  the  free-space  case.  However  the  skin  current  is  still  a  factor  of 
3  lower  than  the  current  caused  by  the  initial  pulse  at  point  3  and  as  a 
result,  we  cannot  consider  that  an  overtest  has  occurred.  The  behavior  of 
the  skin  current  at  late  time  is  very  different  from  the  late  time  behavior 
of  the  skin  current  in  the  free-space  case  indicating  that  the  simulation 
fidelity  is  poor  at  point  10. 

Figure  10  shows  the  normal  component  of  the  electric  field  at 
observer  point  10.  Dramatic  differences  occur  in  the  normal  field  for  the 
strawman  case.  A  large  dip  begins  at  the  time  when  the  wave  is  reflected 
from  the  walls  of  the  tank.  The  component  reaches  a  magnitude  of  approxi¬ 
mately  seven  times  greater  than  in  the  free-space  case  and  has  opposite  sign. 
Furthermore,  the  overall  shape  of  the  curve  is  completely  different.  The 
initial  rise  is  followed  by  a  very  deep  minimum.  Comparison  of  Figures  10c 
and  lOd  indicates  that  the  emission  and  absorption  of  electrons  from  the 
damper  membrane  makes  essentially  no  difference  to  the  F._  field  at  observer 
point  10.  However,  Figure  10c  shows  that  the  presence  of  emission  from  the 
tank  walls  causes  an  enhancement  of  the  reflected  signal.  Evidence  for  this 
same  effect  can  be  seen  in  the  skin  current  plots  (Figure  9)  although  the 
effect  here  is  not  as  pronounced  as  it  is  for  the  electric  field. 

the  reason  lor  this  unexpected  behavior  is  not  fully  understood 
but  a  possible  mechanism  is  suggested  in  Section  1.  The  magnitude  of  the 
fields  at  point  10  is  considerably  less  than  the  magnitude  of  the  fields  in 
the  initial  pulse  at  point  3  so  that  we  cannot  consider  that  an  overtest  has 
taken  place.  The  fields  differ  from  the  frce-spacc  case  to  such  an  extent 
that  the  simulation  quality  criterion  cannot  be  satisfied. 
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3.5 


SATELLITE  OFF  CENTER 


A  scries  of  calculations  was  performed  in  which  the  satellite  was 
removed  from  the  center  of  the  tank  and  placed  close  to  the  rear  wall.  In 
this  configuration,  the  satellite  was  within  36  cm  of  the  damper  membrane. 

The  purpose  of  these  calculations  was  to  examine  the  effects  of  having  the 
satellite  close  to  the  tank  walls.  Stronger  reflections  are  expected  to 
occur  in  this  case  raising  the  possibility  of  overtesting.  The  satellite  is 
now  much  closer  to  the  source  of  electrons  so  that  larger  electric  fields 
on  the  satellite  may  be  produced.  Considerations  of  the  effects  of  the 
proximity  of  the  satellite  to  the  tank  walls  are  of  importance  since  it  is 
probable  that  the  SX'l'F  will  be  used  for  tests  of  very  large  satellites  whose 
dimensions  will  be  comparable  to  those  of  the  tank. 

In  order  to  examine  in  isolation  the  effects  of  position  upon  the 
reflection  of  electromagnetic  waves,  the  fluence  was  not  reduced  from  its 
value  in  the  previous  calculations.  This  assumption  is  not  physically 
justified  since  the  exploding  wire  radiation  behaves  as  a  spherical  wave  falling 
off  as  1/r". 

1 

The  time  integral  of  B“(t)  at  observer  point  3  is  shown  in 
figure  11.  At  observer  point  3,  the  behavior  of  the  integral  is  very 
similar  to  the  behavior  in  the  case  of  the  centered  satellite  (Figure  S). 
figure  11  indicates  that  in  the  case  of  the  shifted  satellite,  the  reflection 
from  the  tank  walls  deposits  more  energy  on  the  satellite  than  it  does  in 
the  case  of  the  centered  satellite  but  the  difference  is  small  (cf.  Figures 
3c  and  11c).  In  the  shifted  satellite  strawman  case,  the  integral  of  B~(t] 
has  a  maximum  value  of  only  10  percent  higher  than  that  attained  in  the 
centered  satellite  strawman  case.  From  figure  lid,  wo  sec  that  the  emitting 
membrane  has  little  effect  on  the  overall  behavior  of  the  integral.  At 
late  time  the  value  of  the  integral  in  the  emitting  membrane  case  is  ~  10 
percent  lowar  than  in  the  M rawmun  case  (compare  Figures  11c  and  lid).  This 
effect  is  similar  t<>  that  -een  in  the  centered  satellite  case  ( compare  Figures 
3c  and  .id  )  . 
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Wall  Emission 


d.  Emitting  and  Absorbing  Membrane 


From  the  plots  in  Figure  11,  we  see  that  the  strawman  case  differs 
from  the  free-space  case  by  less  than  a  factor  of  2  in  /B“(t)dt.  1'his  close 
agreement  between  even  the  most  extreme  cases  indicates  that  in  the  shifted 
satellite  case,  the  simulation  fidelity  at  observer  point  3  is  very  good. 
Further  evidence  for  the  high  quality  of  simulation  is  provided  by  the  plots 
of  skin  current  and  normal  component  of  the  electric  field  shown  in  Figures 
12  and  13  respectively. 

The  general  behavior  of  the  skin  current  is  the  same  in  all  four 
cases;  the  sharp  initial  spike  is  followed  by  lower  amplitude  oscillations. 

In  those  cases  in  which  the  satellite  is  enclosed  in  the  tank,  the  reflected 
pulse  is  quite  prominent  (Figures  12b,  c  and  d)  especially  in  those  cases  in 
which  large  quantities  of  electrons  are  present.  In  no  case  does  the  reflected 
pulse  exceed  ~  80  percent  of  the  initial  pulse  so  that  the  skin  current  plots 
do  not  show  any  evidence  for  the  possibility  of  overtesting.  In  addition, 
Figure  12  indicates  that  the  reflected  pulse  never  provides  skin  currents 
larger  than  twice  the  free  space  values  (compare  Figures  12a  and  12c)  so 
that  the  criteria  for  simulation  fidelity  are  satisfied. 

Comparison  of  the  various  cases  shown  in  Figure  13  show  that  the 
behavior  of  the  electric  field  is  almost  identical  in  all  four  cases.  In 
summary  we  conclude  that  the  influence  of  the  tank  walls  do  not  give  rise 
to  overtesting  and  that  the  quality  of  the  simulations  at  observer  point  3 
is  very  good. 

The  simulation  fidelity  is  much  worse  at  observer  point  10.  The 
plots  of  /B~dt  for  observer  point  10  reveal  that  the  strawman  case  differs 
from  the  free-space  case  by  almost  a  factor  of  200  (Figure  14).  The  presence 
of  the  tank  in  this  case  causes  an  enormous  enhancement  in  the  loule  heating 
of  the  circuitry  at  this  point.  However,  the  maximum  value  reached  by  the 
integral  is  still  comparable  to  the  value  reached  at  point  3  in  the  free- 
space  case  so  that  although  the  simulation  fidelity  is  very  poor,  no  over¬ 
test  i ng  occurs . 


There  are  large  differences  in  the  behavior  of  the  skin  current 
at  observer  point  10.  The  behavior  of  the  skin  current  in  the  strawman  and 
emitting  membrane  cases  (figures  15c,  and  15d)  differs  from  that  of  the 
free- space  case  in  the  same  general  way  as  they  did  in  the  centered  satellite 
cases  (Figures  9c  and  9d)  when  the  satellite  is  close  to  the  tank  walls, 
however,  these  differences  are  much  more  severe.  The  dip  which  occurs  at  the 
time  of  reflection  ~  90  ns  is  approximately  an  order  of  magnitude  larger 
than  the  initial  spike  and  has  an  opposite  sign.  These  large  differences  in 
the  behavior  of  the  skin  current  are  another  indication  of  the  poor  quality 
of  simulation  at  this  point. 

The  reflected  pulse  in  the  perfect  suppression  case  is  much  weaker 
in  the  shifted  satellite  case  than  it  was  in  the  centered  satellite  perfect 
suppression  case  (compare  figures  9b  and  15b.  This  difference  may  be 
attributed  to  the  proximity  of  the  observer  point  to  the  damper  membrane. 

The  most  surprising  results  of  this  work  are  shown  in  Figure  16.  Here 
we  show  the  normal  component  of  the  electric  field  at  observer  point  10. 
reaches  a  magnitude  of  20  times  its  free-space  value  and  has  an  opposite  sign 
in  the  strawman  simulator.  The  magnitude  of  the  field  is  comparable  to  that  at 
observer  3  (sec  figure  13).  The  behavior  seen  here  is  similar  to  but  more  extreme 
than  that  seen  in  the  centered  satellite  case,  with  the  field  reaching  a 
value  of  almost  three  times  greater  than  that  shown  in  Figure  10.  After 
the  dip,  the  field  increases  reaching  positive  values  of  over  six  times  the 
free  space  maximum.  At  late  time  t  >  150  ns,  the  strawman  case  behaves 
differently  from  either  the  free  space  or  perfect  suppression  cases.  Rather 
than  oscillating  about  some  positive  mean  (Figures  16a  and  16b),  the  field 
decreases  monoton  ica  I  I  y  to  very  large  negative  values  (-0.17  esu) .  IVc 
suggest  that  this  decrease  is  caused  by  the  deposition  on  the  satellite  of 
charge  emitted  from  the  tank  walls.  No  significant  differences  were  found 
in  the  behavior  of  the  electric  field  at  observer  point  10  between  the 
strawman  case  and  the  emitting  membrane  case. 
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The  surprising  behavior  of  the  fields  at  observer  point  10  is  not 
well  understood  and  the  next  section  of  this  report  is  largely  concerned 
with  our  attempts  to  understand  this  phenomenon. 


finally,  we  may  draw  the  following  conclusions. 

Simulation  fidelity  is  generally  very  good  at  observer 
points  close  to  the  emitting  surface  of  the  satellite  both 
for  the  perfect  suppression  case  and  for  the  strawman  case. 
Fidelity  is  slightly  worse  when  the  satellite  is  positioned 
close  to  the  tank  walls  but  in  all  cases,  the  criteria  for 
acceptable  simulation  fidelity  are  satisfied. 

In  no  case  was  evidence  for  overtesting  found.  Although 
the  fields  produced  at  observer  point  10  are  much  higher  in 
the  strawman  case  than  they  were  in  the  perfect  suppression 
case,  they  are  never  much  greater  than  they  are  at  observer 
point  3. 

Since  a  satellite  must  be  designed  to  survive  irradiation 
from  any  direction  we  cannot  consider  this  to  be  an  overtest 
although  the  simulation  fidelity  is  clearly  unacceptable. 
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SECTION  4 

ADDITIONAL  CALCULATIONS 


The  peculiar  behavior  of  the  electric  and  magnetic  fields  described 
in  the  previous  section  is  difficult  to  understand.  Because  of  its  implica¬ 
tions  for  the  SXTf  design,  we  feel  that  it  is  important  to  be  confident  of 
these  results  and  to  develop  some  physical  understanding.  The  calculations 
reported  in  this  section  were  aimed  primarily  at  determining  the  cause  of 
the  dip  in  the  magnetic  field  shown  in  figure  15  and  in  the  electric  field 
shown  in  figure  16.  Such  factors  as  the  amount  of  wall  emission  present, 
the  X-ray  spectrum,  the  X-ray  emission  time  history  and  the  method  of  calcu¬ 
lating  the  movement  of  particles  are  examined.  One  purpose  of  these  studies 
was  to  eliminate  the  possibility  that  the  peculiar  behavior  of  the  fields 
is  some  artifact  of  the  code.  The  details  of  the  several  cases  examined  are 
summarized  in  Table  2. 

In  Case  A,  the  satellite  is  in  the  same  position  inside  the  tank 
as  in  the  shifted  satellite  cases  discussed  in  the  previous  section.  Case 
A  differs  from  the  shifted  satellite  strawman  case  only  in  that  electron 
emission  occurs  only  from  the  front  hemisphere  of  the  tank  and  the  damping 
membrane  was  removed.  A  significant  difference  in  the  results  obtained  in 
this  case  is  that  the  normal  component  of  the  electric  field  at  observer 
point  10  has  decreased  by  approximately  a  factor  of  2  (see  figure  17). 

The  dip  seen  in  the  skin  current  in  this  case  (figure  1 H )  is  also  smaller 
than  that  seen  in  the  strawman  case  (figure  15c).  In  addition,  the 
oscillations  in  the  fields  do  not  damp  out  even  at  times  as  late  as  600  ns. 
This  lack  of  damping  is  due  to  the  absence  of  the  damper  membrane.  The 
reduction  in  the  amplitude  of  the  d i p  has  two  possible  causes.  Since  the 
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EW  -  Exploding  Wire 
Br  -  Bremsstrahlung 


Skin  Current  (amps/m)  (statvolts/cm) 


Figure  17.  Case  A— electric  field  at  point  10 


Time  (ns) 

Figure  18.  Case  A— skin  current  at  point  10 
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differences  between  Case  A  and  the  strawman  case  are  in  the  quantity 
and  position  of  the  emitted  charge,  we  conclude  that  the  reduction  in  the 
magnitude  of  the  dip  is  attributable  to  one  or  other  of  these  differences 
or  to  the  combination  of  the  two.  This  calculation  suggests  that  the  dip 
has  its  origin  in  the  sphere  rather  than  in  the  satellite. 

We  next  examine  the  possibility  that  the  peculiar  behavior  of  the 
fields  may  in  some  way  be  caused  by  the  presence  of  the  satellite.  Perhaps 
it  is  some  artifact  of  the  numerical  algorithm  of  treating  the  thin  current 
carrying  struts.  In  order  to  rule  out  this  type  of  possibility,  we  calculated 
Case  B  in  which  the  satellite  was  removed  from  the  tank  and  the  sphere  emission 
is  only  caused  by  the  exploding  wire  X-ray  source.  The  general  behavior  of  the 
electric  field  in  this  case  (shown  in  Figure  19)  is  similar  to  that  found 

for  Case  A  and  to  the  shifted  satellite  strawman  case  shown  in  Figure  16c. 

The  behavior  of  the  skin  current  in  Cases  A  and  B  are  similar  (Figure  20) 
with  roughly  the  same  magnitude  of  dip.  The  reflected  magnetic  field  is 
approximately  40  percent  stronger  in  Case  B  than  in  Case  A  where  its  maximum 
is  0.9  amp/m.  The  initial  spike  seen  at  ~  60  ns  in  Figures  16c  and  17  is 
caused  by  the  interaction  of  the  X-ray  pulse  with  the  satellite.  Since  no 
satellite  is  present  in  Case  B,  no  initial  spike  is  seen.  The  electric 
field  at  observer  point  10  in  Case  B  first  departs  from  zero  at  ~  70  ns  due 
to  the  interaction  of  the  X-ray  pulse  with  the  wall  of  the  tank.  A  clear 
spike  is  produced  reaching  its  greatest  negative  value  at  ~  95  ns.  The 

amplitude  of  this  spike  is  reduced  by  a  factor  of  four  from  the  strawman 

case.  The  absence  of  the  satellite  reduced  the  amplitude  of  the  spike 
but  did  not  cause  it  to  disappear.  Although  the  presence  of 
the  satellite  affects  the  size  of  the  spike,  it  docs  not  cause  it.  We  can 
therefore  exclude  the  possibility  that  the  unexpected  behavior  of  the  fields 
is  some  artifact  of  our  code  description  of  the  satellite.  In  addition,  Case 
B  clearly  demonstrates  that  the  dip  is  in  some  way  caused  by  the  sphere  rather 
than  the  satellite. 
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In  Case  C  we  irradiate  the  sphere  using  the  hremsstrahlung  source 
only.  The  general  behavior  of  the  field  (Figures  21  and  22)  is  remarkably 
similar  to  that  seen  in  the  strawman  case  and  in  Cases  A  and  B  described 
above.  The  magnitude  of  the  dip  is  lower  than  in  Case  A  by  a  factor  of 
two  and  lower  than  the  strawman  case  by  a  factor  of  five.  We  see  that  the 
dip  is  not  strongly  affected  by  the  characteristics  of  the  radiation  striking 
the  sphere. 


In  Cases  I)  and  K,  we  study  the  effect  of  retardation  of  the  pulse. 
The  radiation  source  is  an  exploding  wire  and  particle  emission  is  calculated 
using  a  dipole  driver.  The  dipole  driver  is  placed  on  the  surface  of  the 
sphere  and  has  the  same  magnitude  and  time  history  as  that  placed  on  the 
emitting  surface  of  the  satellite  in  Case  A. 

Case  1)  resembles  Case  A  in  that  we  consider  electron  emission  from 
the  front  half  of  the  sphere  and  the  satellite  is  present  in  the  tank.  Here, 
however,  the  emission  from  the  satellite  is  neglected.  The  resulting  curves 
are  smoother  due  to  the  absence  of  noise  caused  by  particle  statistics. 

Figure  23  shows  the  behavior  of  the  electric  field  in  Case  0.  Because  the 
dipole  driver  on  the  sphere  has  a  large  amplitude  (appropr i ate  in  fact  to 
the  emitting  surface  of  the  satellite),  a  very  strong  dip  occurs  at  the 
time  of  reflection.  The  d i p  is  a  factor  of  two  larger  than  it  is  in  the 
strawman  case.  Following  the  dip,  the  field  rises  to  a  very  large 
positive  peak  (0.03  esu).  The  magnetic  field  at  observer  point  10  (Figure 
211  initially  behaves  very  much  like  the  magnetic  field  in  previous  cases 
(e.g.,  A  and  B)  but  its  behavior  at  late  time  is  more  extreme  reaching 
a  large  positive  peak  and  rapidly  falling  to  a  very  large  negative  value. 

From  case  I'  we  find  that  the  presence  of  a  non-emitting  satellite  has 
little  effect  on  the  dip.  In  addition,  we  arc  able  to  rule  out  the  pos¬ 
sibility  that  the  dip  is  some  artifact  of  our  particle  emitting  and  moving 
code . 
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Time  (ns) 

Figure  22.  Case  C— skin  current  at  point  10 


Figure  23.  Case  D— electric  field  at  point  10 


Figure  24.  Case  D— skin  current  at  point  10. 
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In  Case  I:,  the  satellite  is  absent  and  the  emission  from  the  sphere 
is  from  the  exploding  wire  source  only.  Again  we  use  the  I’  description  of 
the  electron  emission.  In  this  case,  the  dipole  driver  is  not  retarded. 

Instead,  the  dipole  layer  forms  simultaneously  over  the  entire  hemisphere. 
Comparison  of  Figure  25  and  Figure  23  reveals  that  the  dip  in  Case  !•  is 
much  less  pronounced  than  in  Case  1).  The  unretarded  dipole  driver  in  Case  F. 
appears  to  have  considerable  influence  on  the  behavior  of  the  field.  A  local 
minimum  occurs  in  at  the  time  of  reflection,  but  this  seems  to  be  simply 
the  first  in  a  series  of  oscillations.  The  magnetic  field  behaves  similarly 
(Figure  26). 

The  phenomenon  of  the  deep  minimum  is  not  completely  understood. 

IVe  have,  however,  ruled  out  several  possibilities.  This  behavior  has  been 
shown  not  to  be  an  artifact  of  our  code  description  of  the  satellite,  not  an 
artifact  of  our  particle  emission  algorithm  or  of  our  particle  moving 
routines  and  is  not  caused  (although  it  may  be  modified)  by  the  presence  of 
the  satellite.  The  dip  is  not  caused  solely  by  either  the  exploding  or 
hremsst rah  lung  emission.  The  dip  appears  to  he  directly  related  to  the 
emission  area  on  the  sphere  and  is  strongly  influenced  by  the  retardation 
of  the  emission  of  the  sphere. 

One  possibility  is  that  the  dip  arises  from  the  excitation  of  a 
mode  of  the  tank.  The  width  of  the  dip  in  Case  A  is  about  30  ns  suggesting 
that  a  mode  of  frequency  about  33  MHz  may  be  excited.  Iv'e  wish  to  determine 
whether  a  mode  of  such  a  frequency  can  be  excited  by  the  dipole  driver.  In 
order  to  examine  quantitatively  the  excitation  of  the  sphere,  we  have 
calculated  Case  F.  In  this  case  we  remove  the  satellite  and  excite  the 
sphere  using  a  truncated,  unretarded  driver.  The  dipole  driver  is  the  same 
as  that  used  previously  during  the  initial  13  ns  and  has  zero  magnitude  thereafter. 
After  13  ns  the  undriven  sphere  rings  and  its  modal  response  is  shown  on  the 
Fourier  spectrum  of  the  magnetic  field  at  the  center  of  the  mesh  (Figure  27). 
Several  harmonics  are  simultaneously  excited.  One  of  the  strongest  has  a 
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frequency  of  33  MHz  suggesting  that  the  spike  may  be  the  mode  of  the  system 
having  this  frequency.  The  evidence  for  this  is  incomplete  and  many 
difficulties  are  encountered  in  trying  to  explain  the  behavior  of  the  fields 
as  excitations  of  modes  of  the  sphere. 

An  alternative,  perhaps  more  plausible,  hypothesis  for  the  cause 
of  the  dip  will  be  detailed  below. 

Before  describing  the  proposed  hypothesis  for  the  effect  we  will 
try  to  give  some  background  for  understanding  the  behavior  of  the  EM  fields 
in  this  system.  The  effect  of  a  tank  on  an  SGEMP  simulation  can  be  conveniently 
described  in  three  successive  parts.  Prompt  electromagnetic  radiation  effects 
caused  by  electron  emission  from  a  damper  or  tank  walls  arriving  at  the  satel¬ 
lite  at  time,  t,  roughly  equal  to  2R/c,  where  R  is  the  tank  radius  and  c 
is  the  velocity  of  light  (this  time  is  100  ns  for  the  tank  used  in  our 
simulation).  Electromagnetic  radiation  effects  caused  by  reflection  from 
the  tank  walls  can  be  described  in  terms  of  a  modal  picture.  A  modal  descrip¬ 
tion  is  most  convenient  after  the  radiation  has  assumed  a  global  character 
which  occurs  after  t  >  3R/c  when  one  side  of  the  tank  has  felt  the  presence 
of  the  other  side.  Before  this  time  too  many  modes  are  necessary  for  an 
adequate  description.  The  modal  description  is  the  second  part  of  the 
system's  behavior.  The  third  part  arises  from  electrons  emitted  by  the 
walls,  damper  and  source  shields  striking  the  satellite.  Since  electrons 
move  slowly  compared  to  the  velocity  of  light  and  since  the  wall  emission  is 
to  a  first  approximation  spherically  symmetric,  the  effect  of  electron  impact 
can  be  described  by  electrostatics  and  occurs  relatively  late  in  time. 

We  suggest  that  the  spike  or  dip  effect  is  due  to  prompt  electro¬ 
magnetic  radiation.  Recurrence  of  the  effect,  should  modes  add  up  with  the 
proper  phasing,  will  be  diminished  due  to  the  presence  of  the  damper. 
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When  photons  arrive  at  the  sphere  walls  electrons  are  emitted  and 
slowed  by  the  resulting  electric  fields.  These  circumstances  give  rise  to 
hoth  change  of  dipole  moments  and  time  rates  of  change  of  dipole  moments 
at  the  wall.  The  magnetic  field  radiation  from  a  dipole  moment  changing 
in  time  is  given  by 


c 


r  x  p 


(E  =  B  x  r) 


where  r  is  the  vector  from  the  source  to  the  field  point  and  P  is  in  the 
radial  direction.  For  this  geometry  the  prompt  magnetic  field  will  be  zero 
at  the  center  of  the  tank  and  maximum  near  the  edge  where  the  effect  of 
all  the  wall  dipoles  are  additive.  Therefore,  the  hypothesis  compares 
qualitatively  with  the  code  calculation  because  the  dip  effect  is  greatest 
when  the  satellite  is  in  the  shifted  position  and  decreases  by  a  factor  of 
two  when  only  half  the  tank  surface  is  illuminated  or  only  half  the  wall 
dipoles  exist. 


Radiation  from  the  dipoles  do  not  arrive  at  the  satellite 
simultaneously  from  all  wall  positions  so  the  effect  is  spread  out  in  time, 
lor  plane  wave  excitation  the  maximum  time  difference  for  waves  arriving  at 
the  end  of  the  spherical  tank  diametrically  opposite  the  center  of  the  source 
of  the  wave  is  0. 5  R/c.  If  all  points  on  the  wall  are  excited  simultaneously 
the  maximum  time  difference  of  arrival  is  2R/c.  We  should  therefore  expect 
the  ratio  of  the  plane  wave  excited  spike  to  the  simultaneous  wall  emission 
to  be  in  the  ratio  of  4  to  1 .  This  seems  again  to  be  supported  by  the  code 
computations.  A  characteristic  magnitude  for  P  from  the  2  kV  blackbody 
source  is  about  2  x  10^’  statamp/sec  at  the  walls  and  lasting  about  5  ns. 

We  would  therefore  expect  a  prompt  magnetic  field  arriving  at  2R/c  (100  ns) 
of  duration  about  .5  R/c  +  5  ns  =  30  ns  and  having  a  magnitude  of  about 

-7  -y  f  4TrR"  (yyyy)  =  -°35  gauss  ~  3  amp/m  . 
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The  1/2  in  the  above  calculation  arises  because  the  average  of  r  *  I’  is 
estimated  as  1/2;  the  factor  of  5/30  arises  because  the  1’  pulse  is  about 
5  ns  but  is  spread  over  30  ns.  The  estimated  character  of  the  prompt  pulse 
agrees  with  the  observed  character  in  figure  15  and  supports  the  proposed 
hypothesis  for  the  dip  effect. 
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SECTION  5 

RECOMMENDATIONS 


The  presence  of  the  tank  and  electron  emission  from  the  walls  do 
not  appear  to  affect  the  fields  on  and  near  the  satellite  in  such  a  way  as 
to  cause  severe  overtesting.  liven  in  our  most  extreme  case  figure  16c,  the 
magnitude  of  the  field  produced  at  observer  10  is  approximately  the  same  as 
that  seen  at  observer  3.  Our  results  indicate,  however,  that  while 
acceptable  simulation  is  produced  in  the  perfect  suppression  case,  the  more 
realistic  case  of  an  emitting  tank  wall  produces  very  poor  simulation 
fidelity  at  specific  points.  The  departure  of  the  response  from  the  frec- 
space  response  is  most  severe  when  the  satellite  is  placed  close  to  the 
tank  walls,  indicating  that  for  a  strawman  design  SXTF,  serious  fidelity 
problems  may  arise  in  tests  of  larger  satellites. 

Suppression  of  electron  emission  from  the  tank  walls  appear  to  be 
required  to  remedy  this  situation.  Coating  the  tank  walls  with  carbon  (a 
low-photoelectric  yield  substance)  as  in  the  strawman  design  is  inadequate. 
Additional  suppression  or  collimation  is  necessary.  If  collimation  is  not 
sufficient,  an  active  suppression  grid  maintained  at  a  fixed  potential  may 
be  necessary  in  order  to  prevent  wall  electrons  from  causing  a  large  tank 
response . 
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APPENDIX  A 

TREATMENT  OF  THIN  CURRENT  CARRYING  STRUTS 


Our  calculations  attempt  to  describe  electromagnetic  phenomena 
inside  a  tank  of  radius  30  m.  The  test  body  consists  of  three  cylinders 
connected  by  struts  of  radius  1 . 25  cm.  Clearly  it  is  not  possible  to  use  a 
mesh  fine  enough  to  describe  the  struts  while  at  the  same  time  covering  the 
entire  space  of  interest.  Such  a  description  of  the  problem  would  require 
enormous  amounts  of  computer  storage  and  CP  time.  An  alternative  procedure 
is  required. 

Our  technique  of  treating  thin  struts  is  similar  to  that  described 
by  Parks,  et  al.7  and  Lee.8  We  shall  present  a  brief  outline  of  this 
method . 


In  general,  Faraday's  equation  may  be  written 


Integrating  the  equation  and  making  use  of  Stoke's  theorem,  we  obtain 

l/f  • 

Now  consider  this  equation  as  it  applies  to  some  cell  lying  on  the  Z  axis 
and  consider  a  thin  conducting  strut  of  radius  a  <<  Ar  lying  along  the  Z 
axis.  Faraday's  equation  may  then  be  integrated  with  the  area  of  integra¬ 
tion  being  the  shaded  area  in  Figure  A- 1  and  the  contour  being  its  perimeter. 
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Using  the  fact  that  the  field  in  the  vicinity  of  the  strut  varies 
like  1/r,  the  surface  integral  becomes 
Ar/2 


f/^(f)d- 


Carrying  out  the  integration  we  have 


AZ  3B (a) 
c  3t 


a£n(Ar/2a)  . 


The  line  integral  becomes,  since  vanishes  on  the  conductor  and 
the  normal  field  to  a  cylindrical  conductor  varies  as  1/r 


Ar/2 


Ar/2 


-  Fi,  A2  +  E  fa) 


/tdr  -  c(a)r,i /fdr 


Rearranging  the  terms,  the  integrated  Faraday  equation  becomes 


3B(aJ_  1 

3t  c  aHn (Ar/2a) 


r-(a)r,£+l  -  E^rt>. 

AZ 


This  equation  differs  from  the  usual  difference  formulation  of  Faraday's 

equation  onlv  in  that  the  term  in  1:.„  has  been  added. 

L\ 

Ampere's  equation  may  be  treated  similarly.  In  general,  we  have 
1  9f  o  ,  i 

-  rr  =  V  x  ||  -  4ir  ,J  , 

C  a  t  ^  ^ 


where  .1  is  in  abampercs  per  cm“.  Integrating  and  making  use  of  Stoke's 
theorem 


■  da  =  j)  U  ■  dH  -  4tt  f  .J  •  da  , 
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Figure  A-2.  Contour  of  integration  for  H. 


where  the  area  of  integration  is  the  circle  of  radius  Ar  (see  Figure  A-2J 
and  the  contour  of  integration  is  the  circumference  of  the  circle.  Then 
neglecting  the  small  area  of  the  wire 

->  9I-- 

— - =  2irArH,  -  2iTaH(a)  , 

C  dt  1 

which  becomes,  on  rearranging 


3F_ 

“1 

3t" 


which  is  simply  the  usual  difference  formulation  of  Ampere's  equation  with 
an  additional  term. 


The  radial  component  of  Ampere's  law  is 
1  3l;r 

—  =17  x  )1  -  4"JJ  . 
c  .» t  r 


Integrating  over  the  surface  of  the  wire  furnishes  an  equation 
relating  the  normal  electric  field  and  the  surface  magnetic  field. 

These  modifications  of  Maxwell's  equations  permit  the  calculation 
of  the  fields  in  the  vicinity  of  the  wire  and  makes  possible  an  accurate 
description  of  the  electromagnetic  behavior  of  the  test  body  using  a  mesh 
large  compared  to  the  strut  radius.  One  can  thus  realize  large  savings  in 
computer  storage  and  Cl’  time. 
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